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Abstract 

We investigate the effects of the littlest Higgs model with T parity up to the QCD next-to- 
leading order (NLO) on the W^Zjj productions at the CERN Large Hadron Collider (LHC), and 
discuss the kinematic distributions of final decay products and the theoretical dependence of the 
cross section on the factorization/renormalization scale. We find the QCD NLO corrections reduce 
the scale uncertainty of the leading order cross section in case of /if = (J-r- By adopting the 
PROSPINO subtraction scheme (scheme (II)) in analysing the QCD NLO contributions, we can 
obtain the numerical results which keep the convergence of the perturbative QCD description. Our 
results by adopting scheme (II) at the 14 TeV (8 TeV) LHC show that the if -factor for the W£Z H 
production varies in the range of 1.01 ~ 1.10 (1.00 ~ 1.08), while the if-factor for the W^Zh 
production varies in the range of 1.11 ~ 1.13 (1.11 ~ 1-12), when the global symmetry breaking 
scale / goes from 400 GeV to 1.5 TeV (1 TeV). 



PACS: 12.38.Bx, 12.60.Cn, 14.70.Pw 
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I. Introduction 



To interpret the mechanism of electroweak symmetry breaking and resolve the little hierarchy problem 
[T] are the major motivations for the little Higgs models [2J. In those models some new gauge bosons, 
scalars and fermions are introduced at a global symmetry breaking scale / to cancel the one-loop 
quadratic divergences for the Higgs mass from the standard model (SM) O 2] particles. It deserves 
much attention due to their elegant solution to the hierarchy problem and they are proposed as 
one kind of electroweak symmetry breaking models without fine-tuning. Among the little Higgs 
models there is one simplest version, the littlest Higgs (LH) model, providing a set of new heavy 
gauge bosons (Wh,Zh,Ah) and a vector-like quark (T) to implement the divergence cancellation. 
Nevertheless, precision electroweak measurements [5] severely constrain the LH model, especially the 
recent experimental measurements [[51 [7J on the searching for Wjj and Zh bosons. 

The precision electroweak constraints require the LH model characterize a large value of /. To avoid 
fine-tuning between the global symmetry breaking scale / and the electroweak symmetry breaking 
scale, a discrete symmetry named T parity IBJ-fTD] is imposed. In this way, the heavy gauge bosons 
assigned to be T-odd particles do not directly couple with a pair of SM fermions and all dangerous 
tree-level contributions to the precision measurements are forbidden, therefore, the phenomenological 
constraints are somewhat relaxed. Thus the LH model with T parity (LHT) [5]- [T2] deserves more 
attention. In the LHT, heavy gauge bosons, heavy fermions and heavy leptons acquire masses through 
the breaking of the global symmetry, and there exists an attractive dark matter candidate Ah |13j . 
The global symmetry breaking scale / can be lower than 1 TeV [10] , and the processes Wjf — > P v 
and Zh — > are forbidden due to the T parity conservation, leaving the only T-odd heavy gauge 
boson decay modes Wfi —> AhW^ and Zh —> AhH, where H is the lightest neutral Higgs boson, 
followed by the subsequential leptonic decays of and Higgs boson. As a result, the experimental 
constraints f5J [7] on mw H and mz H can not be applied to the T-odd gauge bosons in the LHT. 
Recently, some QCD NLO phenomenological aspects of the LHT have been analyzed in Refs.|14( 115]. 
The WhZh production at the LHC can be significant in searching for the new gauge bosons due to 
the potential of its copious productions as shown in Refs.[HJJ [T7J, where the WhZh production at the 
LHC is only studied at the leading-order (LO). 
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The purpose of this work is to perform a comprehensive analysis for the processes pp — > W^Zh + 
X at the LHC up to the QCD NLO. In Sec. II a brief review of the related LHT theory is given. In 
Sec. Ill we present the details of the calculations. The numerical results and discussions are provided 
in Sec. IV. Finally we give a short summary. 

II. Related LHT theory 

In order to fix notations used in this paper we briefly review the relevant LHT theory. The details of 
the LHT theory can be found in Refs.[8l QUE)]. 

In the LHT the assumed global symmetry SU(5) is broken down spontaneously to SO(5) at 
some high scale / around 1 TeV [18] . Breaking of SU(5) leads to 14 massless Nambu-Goldstone 
bosons, which transform under the electroweak gauge group, SU{2)l x U(1)y, as a real singlet, a real 
triplet, a complex doublet and a complex triplet. Four of the Nambu-Goldstone bosons are treated 
as longitudinal components of the heavy gauge bosons. The others decompose into a T-even SU (2) 
doublet h, identified as the SM Higgs doublet, and a complex T-odd SU(2) triplet <£. 

The T parity transformations for the gauge sector are defined as the exchange between the gauge 
bosons of the two SU(2) x U(l) groups, i.e., VFf o W 2 and B\ <-> B 2 . The gauge couplings of the 
two gauge groups have to be equal, i.e., g\= gi = V2g and g[ = g' 2 = V2g'- Thus their T-odd and 
T-even combinations can be obtained as 

W% = ±={W? - W%), B H = ±(B 1 -B 2 ), (T-odd), 

WZ = -^=(W? + W%), B L = ±(B l + B 2 ), (T-even). (2.1) 
The mass eigenstates of the gauge sector in the LHT are expressed as 

W% = -p(Wfl-=F*W£), Z H = sme H B H + cose H W^, Ah = cobOhBh-suiOhW&, 

Wf = -=(WlTiWl), Z L = -sin6 w B L + cos 9 w Wl, A L = cos 6 W B L + sin 6 w Wl, (2.2) 
v2 



where 9w is the Weinberg angle, and the mixing angle Oh at the C(v| A/ // 2 ) is expressed as 

sin Oh 



±(W-g 12 ) P 



(2.3) 
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The T-even gauge bosons Al, and W L are identified with the SM gauge bosons, while the four 
new heavy gauge bosons, the T-odd partners of SM gauge bosons, are Ah, Zh and W H with masses 
of [16] 

m AH ~ -Lg'f (l - , m ZH ~ m Wll *gf(l- , (2.4) 

where vsm = 246 GeV. The T parity partner of photon, Ah, is the lightest T-odd particle. Therefore, 
the heavy photon is a candidate of dark matter. The masses of SM gauge bosons can be expressed as 
mw = \gvsM and mz = \\/ Q 2 + g' 2 vsM at the tree-level. 

When the T parity is implemented in the fermion sector of the model, the existence of mirror 
partners for each of the original fermions are required. The T-odd partners of SM up- and down-type 
quarks are denoted as t7_ and D_ , where U- = u_ , c_ , t_ and L>_ = d_ , s_ , fe_ . We can get their 
masses as [IB] 

mr/_ ~ v^K/fl-^^f) , m D _=V2 K f, (2.5) 

where k is the mass coefficient in Lagrangian of the quark sector. The Feynman rules in the LHT 
used in this work are presented in Appendix. 

III. Calculations 

In the LO and QCD NLO calculations we employ the FeynArts 3.4 package [19] to generate Feynman 
diagrams and their corresponding amplitudes. To implement the amplitude calculations we apply 
FormCalc 5.4 programs |20| . The t'Hooft-Feynman gauge and the five-flavor scheme (5FS) are adopted 
in this work. 

III..1 LO cross section 

At the parton level the cross section for the qq' — > W^Zh (qq' = ud,us,cd,cs) subprocess in the 
LHT should be the same as that for the corresponding charge conjugate subprocess qq' —> W^Zh 
{qq' = ud, us, cd, cs) due to the CP-conservation. We present the parton level calculations for the 
related subprocess qq' — > W^Zh in this section. By neglecting the contribution of bottom quark in 
the initial state, the LO contribution to the cross section for the parent process pp — > W^Zh + X 
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d Zh d Zh 

(1) (2) (3) 

Figure 1: The LO Feynman diagrams for the partonic process ud — > W^Zh ■ 



comes from the subprocesses 

Q(Pi) + q'(P2) -> W£(p 3 ) + Z H (pi), (qq = ud,us,cd,cs), (3.1) 

where p% (i = 1, 2, 3, 4) represent the four- momenta of the incoming partons and the outgoing , Zh 
bosons, respectively. The Feynman diagrams for the ud — > W^Zh partonic process are shown in Fig{H 
and the LO Feynman graphs for other relevant partonic processes qq' — > W^Zh (qq' = us,cd,cs) are 
similar with those in FigfTJ 

The expression for the LO cross section for the partonic process qq' — > W^Zh (qq' = ud, us, cd, cs) 
has the form as 

&0qg ' = i^4rbl / ? S \M$\ 2 dQ 2 , (qq' = ud, us, cd, cs), (3.2) 

I "I spin color 

where the factors | and | come from averaging over the spins and colors of the initial partons, 
respectively, p is the three- momentum of one initial parton in center-of-mass system, \fl is the partonic 
center-of-mass system energy and is the amplitude of all the tree-level diagrams for the partonic 

process qq' — > W^Zh- The summation is taken over the spins and colors of all the relevant particles 
in the qq' — > W^Zh subprocess. We perform the integration over the two-body phase space of the 
final particles and Zh- The phase space element d£l 2 is expressed as 

Then the LO total cross section for the parent process pp — > W^Zh + X can be expressed as 

cd,cs i i 

a LO = y2 dx i dx 2 [G g /p 1 (xi,fi F )Gg / /p 2 (x 2 ,HF) + (1 ^ 2)] o° ,(s = x x x 2 s), (3.4) 

, 3 -JO Jo 

qq' =ua,us, 



5 



where Gj/p (j = u, c, d, s) is the parton distribution function (PDF) of proton P, which describes the 
probability in finding a parton j with momentum xpj in proton P, s represents the total colliding 
energy squared in the rest frame of proton-proton system, and hf is the factorization scale. 

III.. 2 QCD NLO corrections 

The QCD NLO corrections to the parent process pp — > W^Zh + X at the LHC can be divided into 
four parts: 

• The QCD one-loop virtual corrections to the partonic processes qq' — > W^Zh', 

• The contributions of the real gluon emission partonic processes qq' — > W^Zu + g; 

• The contributions of the real light-(anti)quark emission partonic processes qg — > W^Zh + q'~, 

• The contributions of the PDF counterterms. 

The dimensional regularization method in D = 4 — 2e dimensions is adopted in this work to isolate 
the ultraviolet (UV) and infrared (IR) singularities in the NLO calculations. 

III. .2.1 One-loop virtual corrections to qq' — > W^Zh partonic process 

Some representative Feynman diagrams for the one-loop virtual corrections to the partonic process 
ud — > W^Zh are presented in Figj2j There exist both UV and IR singularities. The masses and 
wave functions of SM quarks and their T-odd partners should be renormalized to remove the UV 
divergences. The counterterms are defined as 

= + (3.5) 
UP* = + (3.6) 
m q- = m q- + 5m q _ , (3.7) 

where t/)q ' R , if)qi R denote the fields of SM quark and T-odd quark, respectively, and m q _ denotes the 
mass of T-odd quark. The on-shell scheme is applied to renormalize the relevant fields and masses, 
then we obtain 

SZi* = [A UV - A IR ] , (3.8) 

6 



(1) 



Z H d 



(2) 



H d d Zh 

(3) 



d 




(4) 




Figure 2: The representative one- loop Feynman diagrams for the partonic process ud — > W^Zh 
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(3.9) 
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where A 



uv 



je + hi(47r) and A/_r 



je + ln(47r). The one- loop virtual contribution 



is UV finite after performing the renormalization procedure. Nevertheless, there still exist soft and 
collinear IR singularities. By adding the contributions of the real gluon/light-(anti)quark emission 
subprocesses and the counterterms of the PDFs at the NLO, the remaining singularities are canceled 
as we shall see later. 

III. .2. 2 Real gluon/light-(anti)quark emission corrections 

The real gluon emission partonic processes for the W^Zh production can be denoted as 

qipx) +q'{P2) W£(p 3 ) + Z H (p4) + g(ps), (qq' = ud,us,cd,cs). (3.11) 

The real gluon emission subprocess qq' — > W^Z^g contains both soft and collinear IR singularities 
which can be conveniently isolated by adopting the two cutoff phase space slicing (TCPSS) method 



In FigJ3]we show the tree level Feynman diagrams for this subprocess. In performing the calculations 
with the TCPSS method, we should introduce two arbitrary small soft cutoff 5 S and collinear cutoff 
5 C . The phase space of the qq' — > W^Zjjg subprocess can be split into two regions: soft gluon region 
(£■5 < ^ S V§) and hard gluon region (E 5 > ^5 S V§) by the soft cutoff 5 S . The hard gluon region is 
separated as hard collinear (HC) and hard noncollinear (HC) regions by the collinear cutoff 5 C . The 
HC region is the phase space where S15 < 5 c s or S25 < S c s (§ij = (pi +pj) 2 ). Then the cross section 
for the real gluon emission subprocess is written as 

2.H ~s , ±.H ±. S , -HC 1 - HC (0 1 o^ 

(Tg = a g + (Tg = O g + (T g . (3.l2) 

According to the Kinoshita-Lee-Nauenberg (KLN) theorem [22], the soft singularity in the soft part 
&g can be canceled by the soft IR divergence in the virtual corrections, while the hard noncollinear 
cross section part a^ c is IR safe. The virtual corrections cancel part of the collinear singularity and 
the PDF counterterms absorb the remaining collinear divergence. 

Beside the real gluon emission subprocesses, the real light-(anti)quark emission subprocesses, which 
have the same order contributions with the real gluon emission subprocesses, should be taken into 
account. This kind of subprocesses is denoted as 

?(Pl) +g(P2) -> W^(p 3 ) + Z H (p A ) +q'(p 5 ), {qq =ud,cs,du,sc). (3.13) 

The corresponding Feynman diagrams for the subprocess qg — > W^Zh + q' at the tree-level are shown 
in FigHJ Using the TCPSS method described above, the phase space can be split into a collinear 
(C) region (,§15 < 5 c s or S25 < S c s) and a noncollinear (C) region (§15 > 5 c s and §25 > 5 c s) by a 
collinear cutoff S c . Therefore, the cross section for the real light- (anti) quark emission subprocess can 
be expressed as 

a? = a° + af. (3.14) 

The cross section in the noncollinear region is finite and can be evaluated in four dimensions using 
Monte Carlo technique while a^j contains collinear singularity. After adding the renormalized virtual 
corrections and the real gluon/light-(anti)quark emission corrections to the subprocess qq 1 — > W^Zh , 
the partonic cross section still contains the collinear divergence, which can be absorbed into the 
redefinition of the PDFs at the NLO. 
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Figure 3: The tree-level Feynman diagrams for the real gluon emission partonic process qq' 
WftZu + g, (qq' = ud, us, cd, cs). 



III. .2. 3 PDF counterterms 



The PDF counterterms, 6Gi/p(x, hf) (i = 9,u,u,d,d,c,c,s,s), which absorb the remaining collinear 
divergence, can be split into two parts: the collinear gluon emission part 6G^ l p° n \x, /if) and the 
collinear light-quark emission part 5G^j p ark \x, uf)' 



6G q(g)/p( x ^F) = 5G^ } ° p \x, fi F ) + 8G ( ^yp> (x, (j, F ), (q = u,u,d,d,c,c,s,s), (3.15) 



i(quark) 



where 



6G ( °£™ p \x,iif) = - 
5G^ h \x,n F ) = \ 



a s r(l -e) ( 



4 xe 



5G^ k \x^ F ) 



2vrr(l - 2e) V m| 
' a s T(l - e) / 4-Kfi R 
27rr(l-2e) V P? F 

a s T(l -e) f 4vr/i| 



dz 



Pq q (gg)(z)G q (g)/p{x/z,V F ), 



2 \ e-i r l 



dz 

—P qg (z)G g/P (x/z,fJ,F), 



27rr(l - 2e) \ m| 



e-i c,c,s,s,b,b, ^ 



Pgq{z)G q/ p(x/z,fl F ). (3.16) 



More details about the explicit expressions for the splitting functions Pij(z)(ij = qq,qg, gq, gg) are 
available in Ref. 21 . 



III. .2.4 Total QCD NLO correction 

Finally, we have eliminated all the UV and IR singularities by performing the renormalization proce- 
dure and adding all the QCD NLO correction components, and we get the finite QCD NLO corrected 
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Figure 4: The tree-level Feynman diagrams for the real light-quark emission partonic process qg —> 
W^jZh + 4 (qq' = u d, cs, du, sc). 

integrated cross section for the pp — > W^Zh + X process as 

onlo = (7 LO + Aa NLO = (7 L o + Aa (2) + Aa (3) . (3.17) 

The two-body term A<j( 2 ) includes the one-loop corrections to the pp — > W^Zh + X process and 
the tree-level contributions in the soft and hard collinear regions for the real gluon/light-(anti)quark 
emission processes, while the three-body term A<j( 3 ) contains the cross sections for the real gluon/light- 
(anti)quark emission processes over the hard noncollinear region. 

In this work, two event selection schemes are applied in discussing the QCD NLO corrections. 
In scheme (I) all the NLO correction components mentioned above are included in the QCD NLO 
corrections, called also the inclusive event selection scheme. In this scheme, there exists resonance 
effect in FigsJU (5)-(8) duo to the possible on-shell (/_ propagator and those Feynman diagrams could 
lead to large corrections to the Born pp — > W^Zh + X process, so that the perturbative convergence 
would be eventually destroyed. To deal with the resonance effect in these partonic processes, the q- 
mass squared rr? q _ in its propagator should be replaced by rr? q _ — im q _T q _. The partial decay widths 
of T-odd quarks are obtained numerically by adopting the expressions presented in Ref . [15] . 

Actually, the contributions from the diagrams for the qg — > W^Z^ + q 1 subprocess with intermedi- 
ate on-shell T-odd quark g_ shown in FigsJU (5)-(8), should pertain to other on-shell Wnq~ and Znq~ 
associated production channels, i.e., pp — > qg — > Wnq'- + X and pp — > qg — )• Zuq~ + X processes, 
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followed with subsequential decays of q'_ — > Z^q' and g_ — > Wnq 1 , respectively. To avoid double 
counting and to keep the convergence of the perturbative QCD description for the pp — > W^Zh + X 
process, we adopt the PROSPINO subtraction strategy [23, 24J to remove the on-shell T-odd quark 
q- contributions called scheme (II). This subtraction scheme can provide a reliable production rate 
since it only subtracts the squared on-shell amplitudes and does this point by point over the entire 
phase space. The PROSPINO subtraction is done by performing a replacement of the Breit-Wigner 
propagator 

\M\ 2 (sy Hq ) ^ \M\ 2 (sy Hq ) 

{s Vm - m\_) 2 + m\_Y\_ {s VH q - m\_) 2 + m\_Y\_ 

\M\ 2 (m\ \ 



{svm _ m 2 )2 + m 2 W ^ " 4 ^-)©K- " ( 3 - 18 ) 



where sy Hq is the squared momentum flowing through the intermediate g_ propagator. 

Analogously, we can follow above calculation procedure to evaluate the LO and NLO QCD cor- 
rected results for the pp — > W^Zh + X process at the LHC. 

IV. Numerical results and discussions 
IV..1 Input parameters 

The two mixing matrices, Vh u and Vnd cannot be set to be unit matrices simultaneously due to the 
condition of V^ u V Hd = Vckm [25] . In our numerical calculations Vhu 

is set unit matrix, then 

we get V H d = Vckm- We take a cw (m 2 z )' 1 = 127.916, m w = 80.399 GeV, m z = 91.1876 GeV 

= 1 - /W n2 



and sin 0w = 1 — \~^) = 0-2226 [26J. We neglect the masses of /u-lepton and light quarks. The 
colliding energy in the proton-proton center-of-mass system is set as yfs = 8 TeV for the early LHC 
and y/s = 14 TeV for the future LHC. We define (jlq = (mw H + m z H )/2 and adopt CTEQ6L1 and 
CTEQ6M PDFs in the LO and NLO calculations, respectively. The LHT T-odd quark mass coefficient 
parameter k is fixed to be 1. Consequently the masses of heavy gauge bosons and T-odd quarks are 
only the functions of the LHT parameter / as shown in Eqs. (|2.4p and (|2.5p . The Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements are taken as 

0.97418 0.22577 
Vckm = j V cd V cs V cb j j -0.22577 0.97418 | . (4.1) 

1 



v ud 


V us 


v ub 


Vcd 


Vcs 


V cb 


v td 


Vts 


v tb 
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/ 


m WH « m ZH 




m u _ = m c _ 


m d _ = m s _ 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


500 


322.1 


67.5 


685.7 


707.1 


700 


457.8 


102.7 


974.7 


989.9 


800 


525.1 


119.7 


1118.0 


1131.4 


900 


592.3 


136.4 


1260.9 


1272.8 


1000 


659.3 


153.0 


1403.5 


1414.2 


1100 


726.1 


169.4 


1545.9 


1555.6 


1300 


859.7 


202.0 


1830.3 


1838.5 


1500 


993.1 


234.5 


2114.2 


2121.3 



Table 1: The masses of Wh, Zjj, Ah and g_ (q_ = u_, d_, c_, s_) for some typical values 
of the LHT parameter / with k = 1. 

By using Eqs. (|2.4H2.5p and taking the LHT parameter k = 1, we obtain the masses of heavy gauge 
bosons and T-odd quarks for some typical values of the LHT global symmetry breaking scale / and 
list them in Table [TJ 

IV.. 2 Checks 

The correctness of our calculations are verified through the following aspects: 

1. Our LO cross sections are in good agreement with the results read out from Fig. 9 of Ref.|16j 
when we employ the same input parameters and PDFs as used in Ref.|16|. 

2. After combining all the contributions at the QCD NLO, the cancelations of UV and IR diver- 
gences are verified. 

3. We make the verification of the S s /S c independence of the total QCD NLO correction, where 
two arbitrary cutoffs 5 S and 5 C [21] are introduced to separate the phase space in order to isolate the 
soft and collinear IR divergences, respectively. Eq. (|3.17p shows that the total QCD NLO correction 
(Actnlo) is obtained by summing up the two-body and three-body corrections (A<j( 2 ) and Acr^ 3 )). We 
depict A<j( 2 \ Ac - ® and Aa^LO for the process pp — > ud — )• W^Zh + X as functions of the soft cutoff 
5 S in FiglSJa) with / = 600 GeV, k = 1, 8 C = <5 S /100 and fi = fi = (m W[1 + m z J/2 = 390.20 GeV. 
The amplified curve for the total correction Aonlo in Figl5ja) is demonstrated in FigO^b) together 
with calculation errors. From these two figures we find that the total QCD NLO correction A.unlo 
is independent of the two cutoffs within the statistical errors. This independence is an indirect check 
for the correctness of our work. We adopt also the dipole subtraction (DPS) method [27] to deal with 
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the IR singularities. The total QCD NLO correction Aotvlo obtained by adopting the DPS method 
with ±lcr statistic error is plotted as the shadowing region in FigUJb). We can see that the results 
from both the TCPSS method and the DPS method are in good agreement. In further numerical 
calculations, we fix 5, = 1 x 10 -4 and 5 r = 1 x 10 -6 . 





Figure 5: (a) The dependence of the QCD NLO corrections to the pp — > ud — >■ W^Zjj + X process on 
the cutoffs 5 S and 5 C at the yfs = 14 TeV LHC, where we take / = 600 GeV, k = 1, S c = 5 S /W0 and 
A 4 = Mo = 390.20 GeV. (b) The amplified curve for Aatot in FigJS^a). The shadowing region shows 
the result by adopting the DPS method with ±1<t statistic error. 



IV. .3 Dependence on factorization/renormalization scale 

In order to investigate whether the production rates for the pp — > W^Zh + X and pp — > WJjZh + X 
processes at the y/s = 14 TeV and the \fs = 8 TeV LHC have the stabilization of the dependence on the 
unphysical renormalization and the factorization scales, we present FigsJUJ^a) and (b) to describe the 
cross sections as functions of the renormalization and the factorization scales varied independently and 
simultaneously. We show the cross section profile both at the LO and at the QCD NLO by adopting 
the event selection scheme (II) and taking the LHT parameters / = 1 TeV and k = 1. The curves of the 
(Tlo and (Jnlo for the pp — > W^Zh + X process are labeled by "LO+" and "iVLO + ", while those for 
the pp — > WflZn + X process are labeled by "LO_" and "NLO-" , respectively. The two figures trace 
the scale dependence following a contour in the [ir — [ip plane as shown in each left panel of Figs ^ a) 
and (b). From these figures we can see that the QCD NLO corrections do not obviously improve 
the scale uncertainty with individual variation of either [ir or fip. Particularly, the LO partonic 



13 



processes for thepp-> W±Z H + X processes are pure electroweak channels where the (j,r dependence 
is invisible at the LO, as shown in FigsGD(a)-(3), (a)-(5), (b)-(3) and (b)-(5). FigslBJa)-(l) and (b)-(l) 
show that the scale uncertainty is reduced by the NLO corrections with simultaneous variation of 
and fip. It demonstrates that when we set = \ip and vary both scales simultaneously, it may 
lead to artificial cancelations among renormalization and factorization logarithms, and thus hiding 
the scale dependence. In the following discussions the factorization/renormalization scale is fixed as 
IM) = {mw„ +m ZlI )/2. 




2.0 

O 1.5 



(1) 



(2) PR = 5fi 



W+ 



(3) J u F =0.2 rt) 



(4) W =0.2 W) 



(5) [if = 5fi a 

v/i = 8 TeV 
K = 1 

/ = 1 TeV 



Mr,f/po 



Figure 6: Profile of the renormalization and factorization scale dependence of the LO and NLO 
corrected cross sections for the processes pp — > W^Zh + X and pp — > WJjZh + X. The two plots 
trace the scale dependence following a contour in the [ir — \ip plane. There we take the (II) selection 
scheme and assume /i/Vo = 0.2 — 5, the LHT parameters / = 1 TeV and k = 1. (a) at the -^s = 14 TeV 
LHC. (b) at the ^7s = 8 TeV LHC. 



IV.. 4 Dependence on global symmetry breaking scale / 

We depict the LO, QCD NLO corrected integrated cross sections and the corresponding K-factors for 
the pp — > W^Zh + X and pp — > WJjZh + X processes as functions of the global symmetry breaking 
scale / at the y/s = 14 TeV and the y/s = 8 TeV LHC in FigsJ7J[a), (b), (c) and (d), respectively, 
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with k = 1. The curves labeled by "NLO I" and "NLO II" are for the QCD NLO corrected cross 
sections using the (I) and (II) selection schemes, respectively. FigsJ7]^a,b,c,d) demonstrate that the LO 
and QCD NLO corrected total cross sections for the pp — > W^Zh + X processes decrease sensitively 
with the increment of / due to the fact that the masses of final Wjj and Zh become heavier and 
consequently the phase space becomes smaller as the increment of /. The numerical results for the 
pp — > W^Zh + X processes at the LHC for some typical values of / are presented in Table [2J 
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Figure 7: The LO, QCD NLO corrected integrated cross sections and the corresponding K-factors 
as the functions of the global symmetry breaking scale / with k = 1. (a) for the pp — > W^Zh + X 
process at ^ = 14 TeV LHC. (b) for the pp W^Z H + X process at Jl = 8 TeV LHC. (c) for the 
PP -> WJjZ h + X process at ^/s = 14 TeV LHC. (d) for the pp^W^Z H + X process at *fs = 8 TeV 
LHC. 
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(TeV) 


/ 

(GeV) 


a LO 
(fb) 


°NLO 
(fb) 


k( w+ ) 


(fb) 


AW-) 

NLO 


(fb) 




500 


321.096(8) 


350.8(1) 


1.09 


136.130(5) 


153.12(7) 


1.12 




700 


72.055(2) 


77.26(3) 


1.07 


26.888(1) 


30.03(1) 


1.12 


14 


900 


21.9589(5) 


23.159(7) 


1.05 


7.3867(3) 


8.216(3) 


1.11 




1100 


7.8997(2) 


8.203(3) 


1.04 


2.44038(9) 


2.709(1) 


1.11 


8 


500 
700 
900 


94.168(2) 
15.7549(4) 
3.49785(8) 


100.01(5) 
16.259(7) 
3.515(1) 


1.06 
1.03 
1.01 


32.379(1) 
4.6757(2) 
0.93703(3) 


36.10(4) 
5.192(6) 
1.044(1) 


1.11 
1.11 
1.11 



Table 2: The numerical results of a^ Q , cr^^ , ^ lo , ct^ lo for the pp — > W^Z^ + X 
and pp — > WJjZh + X processes and their corresponding K- factors at the -y/i = 14 TeV 
and the ^fs = 8 TeV LHC by adopting the event selection scheme (II) and taking k = 1, 
/i = /io for some typical values of /. 



IV.. 5 Differential cross sections 



In this subsection we focus on the kinematic distributions of final decay products. The WhZh associ- 



ated production at the LHC are followed by the heavy gauge boson decays of — > W t Ah — > /U T 



(-) 



Ah and Zh — > HAh- The branching ratios of decays for the Wh boson, Zh boson and W boson are 
taken as Br(W H -»■ WA H ) = 100%, Br(Z H -»■ HA H ) = 100% for k = 1 and / = 1 TeV p] and 
Br(W T — > ^ = 10.57% |26| . respectively. In the following we consider the WhZh production 



channel including its subsequential decays as 



pp 



(-) 



W%Z H -> W^A H HA H -> f V A H HA H 



(4.2) 



Thus one expects that the WhZh production at the LHC could be detected via the fi T H + $ T ($ T 
= transverse energy of Vu +2Ah) channel. 

The LO, QCD NLO corrected transverse momentum distributions of W boson and the light neutral 
Higgs boson H for the pp — > W^Zh + X and pp — > W^Zh + X processes, and the corresponding 
K-factors in scheme (II) at the yfs = 14 TeV LHC and the y/s = 8 TeV LHC are presented in 
Figs|8ja,b,c,d) and Figsj9ja,b,c,d) separately. There we take / = 1 TeV and n = 1. From these four 
figures we find that the QCD NLO corrections enhance the LO transverse momentum distributions in 
most plotted ranges of pt, and the iT-factors are all less than 1.20. The maxima of the distributions 
drJh °$ LO and daL °^ LO are all located at about p T ~ 220 GeV. 
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Figure 8: The LO, QCD NLO corrected pYp distributions and the corresponding K-factors in scheme 
(II) for the pp ^W^Z H -»• W ± A H HA H + X processes by taking / = 1 TeV and k = 1. (a) for the 
pp -»• WjZn -»• VF+Ah^Ah + X process at the ^/i = 14 TeV LHC. (b) for the pp -»• WjZ H -»• 
VF + Ah#A h + X process at the = 8 TeV LHC. (c) for the pp -> -> W^-A^iTA^ + V 

process at the ^/s = 14 TeV LHC. (d) for the pp -)• W^Z H ->■ W~A H HA H + V process at the 
Vs = 8 TeV LHC. 
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Figure 9: The LO, QCD NLO corrected ptf distributions and the corresponding i^-factors in scheme 
(II) for the pp -> W^Z H -»• W ± A H HA H + X processes by taking / = 1 TeV and k = 1. (a) for 
the -»• VF+Z^ -»• W+Ajyif + X process at = 14 TeV LHC. (b) for the pp -> W£Z H -»• 
W + A H HA H + X process at ^/i = 8 TeV LHC. (c) for the pp -> W H Z H -»• W~A H HA H +X process 
at = 14 TeV LHC. (d) for the pp -»• W^Z H -> W-Ajyif A H + V process at y/s = 8 TeV LHC. 
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The LO and QCD NLO corrected transverse momentum distributions of the final /i-lepton and 
missing energy (AhAh i/^) for the pp — > W^Zh — > h~vAhHAh + X and pp — > W^Zh — > 
[i^v^AhKAh + X processes, and the corresponding ^-factors in scheme (II) at the early LHC and 
the future LHC are depicted in Figs ITOT a.b.c.d) and Figs fTTT a.b.c.d). respectively. There we take 
/ = 1 TeV and k = 1. FigsfTOTa.b) are for the px distributions of /i~, and FigsJTUTc.d) for fi + , 
respectively. FigsJIUI (a), (b), (c) and (d) demonstrate that both the LO and the QCD NLO corrected 
p^ distributions at both the early LHC and the future LHC decrease rapidly with the increment of 
pt^. FigsJll(a.b.c.d) show that the LO and NLO missing transverse momentum distributions reach 
their maxima at p™ tss ~ 290 GeV. 

To show how the 0{a s ) contributions correct the LO differential cross sections at the future and 
early LHC, we depict the LO, QCD NLO corrected rapidity distributions of final VF-boson and Higgs 
boson (\yw\ and \vh\) for the W^Zh production processes in Figs fTTT a.b.c.d) and Figs fT5r a.b.c.d). 
respectively. The W + and Higgs boson rapidity distributions of the pp — > W^Zjj — > W + A^H Ah + X 
process at the future and early LHC are depicted in Figs fT27 a.b) and FigsJTSTa.b). respectively. 
FigsJT2Tc.d) and FigsJTSTc.d) provide the and \yu\ distributions of the pp — > W^Zh — > 

W~ AhHAh + X process, which offer the comparisons with Figsll2(a.b) and Figs]13fa.b) correspond- 
ingly. The rapidity distributions of the final /i-lepton (|y^|) at the LO and QCD NLO are presented 
in Figs mT a.b.c.d). The y^+ distributions for the pp — > W^Zh — > ^u^AhHAh + X process at the 
y/s = 14 TeV and yfs = 8 TeV LHC are plotted in Figs fHT a.b) separately, while the y^- distributions 
for the pp -)■ WjjZ H -)■ iTv^AhHAh + X process at the y/s = 14 TeV and y/s = 8 TeV LHC 
are shown in Figs[HJc,d) respectively. All these figures are obtained by taking the LHT parameters 
/ = 1 TeV, k = 1 and adopting the event selection scheme (II). The corresponding K- factors are also 
plotted in each nether plot of Figsll2|. Figsll3l and Figsll4l We can see from all these figures that the 
QCD NLO corrections do not make shape change in the rapidity distributions. 

V. Summary 

We present the calculations of the WhZh production at the CERN LHC up to the QCD NLO 
in the littlest Higgs model with T parity. The dependence of the cross section on the factoriza- 
tion/renormalization scale are investigated theoretically, and the rapidity and transverse momentum 
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Figure 10: The LO, QCD NLO corrected distributions and the corresponding K- factors in scheme 

(II) for the pp — > W^Zh — > ^ ^ AhHAh + X processes by taking / = 1 TeV and k = 1. (a) for 
the pp ->• W^jZ H ->• [i + v^A H HA H + X process at a/s = 14 TeV LHC. (b) for the pp ->• W^Z H ->• 
/i + i/ M ^4fl-fL4f/+V process at ^/s = 8 TeV LHC. (c) for thepp — > WJjZh — > fi" u^AhH Ah+X process 
at yfs = 14 TeV LHC. (d) for the pp -^W^Z H -> ^D^AhHAh + X process at y/s = 8 TeV LHC. 



20 



5.0 




1.4 F ' ' ' ' ' ' ' ' 1 1.4 F 




300 600 900 300 600 900 

pf ss (GeV) p™' ss (GeV) 



Figure 11: The LO, QCD NLO corrected p™ lss distributions and the corresponding ^-factors in 

scheme (II) for the pp — > W^Zh — > ^ AhHAh + X processes by taking / = 1 TeV and 
k = 1. (a) for the pp ->• ->• ^l + u^A h HA h + X process at the a/s = 14 TeV LHC. (b) 

for the pp ->■ W^Zh — > h+v^AhHAh + X process at the \fs = 8 TeV LHC. (c) for the ->■ 
W^Z H ->• pTv^A H HA H + V process at the ^/i = 14 TeV LHC. (d) for the pp ->• W^Zh ->• 
11~v^AhHAh + V process at the y/s = 8 TeV LHC. 
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Figure 12: The LO, QCD NLO corrected rapidity distributions of the of final IF-boson, \yw\-, an d 
the corresponding if- factors in scheme (II) for the pp — > W^Zh — > W^AhHAh + X processes 
by taking / = 1 TeV and k = 1. (a) for the pp ->■ W^Z H -»■ W + A H HA H + X process at the 
Vs = 14 TeV LHC. (b) for the W^Z H W + AhH Ah + X process at the v/s = 8 TeV 

LHC. (c) for the pp W~A H HA H + X process at the y/s = 14 TeV LHC. (d) for the 

-> W^Zh ->• W~A H HA H + V process at the ^/i = 8 TeV LHC. 
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Figure 13: The LO, QCD NLO corrected corrected rapidity distributions of final Higgs boson, \uh\, 
and the corresponding ET-factors in scheme (II) for the pp — > W^Zh — > W^AhHAh + X processes 
by taking / = 1 TeV and k = 1. (a) for the pp ->■ W^Z H -»■ W + A H HA H + X process at the 
Vs = 14 TeV LHC. (b) for the pp -»• W^Zjf -»• W + AhH Ah + X process at the v/s = 8 TeV 
LHC. (c) for the pp -»• W^Z/r -»• W-AfffTAjj + V process at the y/s = 14 TeV LHC. (d) for the 
pp ->• W# Zjj ->• W~A H HA H + V process at the ^/i = 8 TeV LHC. 
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Figure 14: The LO, QCD NLO corrected rapidity distributions of final muon lepton and the cor- 
responding K-factors in scheme (II) for the pp — > W^Zh — > ^ AhHAh + X processes by 
taking / = 1 TeV and k = 1. (a) for the pp — > W^Zh — > ^v^AfjHAn + X process at the 
v/s = 14 TeV LHC. (b) for the pp W%Z H ^ + v^A H HA H + X process at the y/s = 8 TeV 
LHC. (c) for the pp W H Z H ^-u^A H HA H + X process at the a/s = 14 TeV LHC. (d) for the 
pp -)• WJjZ h -)• yrv^AuHAu + V process at the ^/i = 8 TeV LHC. 
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distributions of final decay products at both LO and NLO are presented. For the purpose of provid- 
ing reliable predictions on the pp — > W^Zh + X process at the LHC, we adopt two event selection 
schemes in considering the QCD NLO corrections for comparison. By using the inclusive scheme 
the perturbative convergence could be destroyed, while we can keep the convergence of the perturba- 
tive QCD description and get moderate QCD NLO corrections to the production rate with evidently 
reduced scale uncertainty by adopting the PROSPINO subtraction scheme and setting \xf = ^r- 
With this scheme the QCD NLO correction enhances the LO cross section, and the corresponding K- 
factor for the W^Zfj production process at the future (early) LHC varies in the range of 1.01 ~ 1.10 
(1.00 ~ 1.08) when / goes from 400 GeV to 1.5 TeV (1 TeV), while the if-factor for the W H Z H 
production process at the future (early) LHC varies in the range of 1.11 ~ 1.13 (1.11 ~ 1.12) in the 
same region of /. 
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VI. Appendix 

We list the Feynman rules for the coupling vertices in the LHT related to this work in Table [3] 
P [TCI [281 E9], where P LjR = ±(1 T 7s) and v = v SM . 



Vertex 


Feynman rule 


Vertex 


Feynman rule 


<5 i _D j (i,j = l,2) 
ZH^Ui-Uj&j = 1,2) 


[9""(fcl - k 2 )" + 

a vp (k 2 - fe 3 ) M + 9 w (fc 3 - fcl)"] 

i-^l^PL(V Hd )ij 


92 ^ g; 

^A-CjliJ = 1.2) 
Z H n£>i-D s (i,i = 1,2) 


!9 s (T a )„ / 37 fI 



Table 3: The related LHT Feynman rules used in this work, Ui = u,c, Di = d,s, Ui- = U—,c- and 
L>j_ = d- , s_. i is the generation index. 
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